Detecting single-electron tunneling involving virtual processes in real time 
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We use time-resolved charge detection techniques to probe virtual tunneling processes in a double 
quantum dot. The process involves an energetically forbidden state separated by an energy 5 from 
the Fermi energy in the leads. The non-zero tunneling probability can be interpreted as cotunneling, 
which occurs as a direct consequence of time-energy uncertainty. For small energy separation the 
electrons in the quantum dots delocalize and form molecular states. In this regime we establish the 
experimental equivalence between cotunneling and sequential tunneling into molecular states for 
electron transport in a double quantum dot. Finally, we investigate inelastic cotunneling processes 
involving excited states of the quantum dots. Using the time-resolved charge detection techniques, 
we are able to extract the shot noise of the current in the cotunneling regime. 



I. INTRODUCTION 

A semiconductor double quantum dot (DQD) is the 
mesoscopic analogue of a diatomic molecule. The energy 
levels and the interdot coupling energy can be precisely 
controlled with gate voltages [l], which allows the DQD 
to be tuned to a configuration where the electron wave- 
functions hybridize and form molecular states extending 
over both QDs. The DQD thus provides a tunable two- 
level system, which has been utilized to perform coherent 
manipulation of a single charge in semiconductor nanos- 
tructures [H [3] . 

An alternative approach to molecular states at large 
detuning is to study electron transport in the DQD in 
the framework of cotunneling [4 . Cotunneling involves 
an electron (or hole) that virtually tunnels through an 
energetically forbidden charge state of the QD positioned 
at an energy 5 away from the Fermi energy in the leads. 
The process occurs on a timescale Tcot ^ T^/^ limited 
by time-energy uncertainty [5 . Cotunneling currents are 
generally small and difficult to measure, but the effect 
has been utilized for QD spectroscopy [H [7] , for studying 
cotunneling- mediated transport in single QDs [8], or for 
investigating spin effects in double QDs [9^ . 

In this work we use a quantum point contact (QPC) 
as a charge sensor [10 to detect single-electron tunneling 
in the DQD in real-time [111 |12] [13^ . Similar setups have 
been used for investigating single-spin dynamics [M] , de- 
tecting single-particle interference [15], probing interac- 
tions between charge carriers [16] or for measuring ex- 
tremely small currents [T7l[T8]. Here, we utilize the tech- 
nique to count electrons cotunneling through the DQD. 
The method provides a precise measurement of the tun- 
neling probability as a function of energy separation 5 
between the QDs, allowing a direct comparison with the 
rate expected from time-energy uncertainty. In the limit 
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of 5 ^ 0, the electrons form molecular states extend- 
ing over both QDs. Here, we measure tunneling rates 
expected from sequential tunneling into bonding and an- 
tibonding states of the DQD. The results experimentally 
establish the equivalence between cotunneling into cou- 
pled QD states and sequential tunneling into molecular 
states of the DQD. 

In the last part of the paper we investigate inelastic co- 
tunneling processes involving excited states of the DQD. 
Finally, we use the time-resolved charge-detection tech- 
niques to extract the shot noise of the DQD current in 
the cotunneling regime. 



II. EXPERIMENTAL SETUP AND METHODS 

The measurements were performed on the structure 
shown in Fig. [ija). The sample is fabricated with lo- 
cal oxidation [19j of a GaAs/Alo.sGao.TAs heterostruc- 
ture, containing a two-dimensional electron gas (2DEG) 
34 nm below the surface. The structure consists of two 
QDs (marked by 1 and 2 in the figure) connected by two 
separate tunnel barriers. Each QD contains about 30 
electrons. For the results presented here only the upper 
tunnel barrier was kept open; the lower was pinched- 
off by applying appropriate voltages to the surrounding 
gates. The sample details are described in Ref. [20 . 

The electron population of the DQD is monitored by 
operating the QPC in the lower-right corner of Fig. [TJa) 
as a charge detector [10]. By tuning the tunneling rates 
of the DQD below the detector bandwidth, charge tran- 
sitions can be detected in real-time [HI |T2l |T3] . In the 
experiment, the tunneling rates Fs and Fd to source and 
drain leads are kept around 1 kHz, while the interdot cou- 
pling t is set much larger {t ~ 20 /ieV ~ 5 GHz). Interdot 
transitions thus occur on timescales much faster than it 
is possible to register with the detector (rdet ^ 50 /is) 
^ , but the coupling energy may still be determined 
from charge localization measurements [22 • The con- 
ductance of the QPC was measured by applying a bias 
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FIG. 1: (a) AFM-image of the sample. The structure con- 
sists of a double quantum dot (DQD) (marked by 1 and 2) 
with a near-by quantum point contact, (b) Charge stability 
diagram of the DQD, measured by counting electrons enter- 
ing and leaving the DQD. The data was taken with a voltage 
bias of Vdqd-sd = 600 /xV applied over the DQD. The QPC 
conductance was measured with Vqpc-sd = 300 /xV. The 
count rates were extracted from traces of length T = 0.5 s. 
(c) Energy level diagrams for different configurations in (b). 



voltage of 200 — 400 /iV and monitoring the current [/qpc 
in Fig.jlja)]. We ensured that the QPC bias voltage was 
kept low enough to avoid charge transitions driven by 
current fluctuations in the QPC [20]. The sample is real- 
ized without metallic gates so that the coupling between 
dots and QPC detectors is not screened by metallic struc- 
tures. 

Figure [ijb) shows a charge stability diagram for the 
DQD, measured by counting electrons tunneling into and 
out of the DQD. The data was taken with a bias voltage 
of 600 /iV applied across the DQD, giving rise to finite- 
bias triangles of sequential transport [1 . The diagrams 
in Fig. [TJc) show schematics of the DQD energy levels 
for different positions in the charge stability diagram. 
Depending on energy level alignment, different kinds of 
electron tunneling are possible. 

At the position marked by I in Fig. [ijb), the electro- 
chemical potential jii of QDl is aligned with the Fermi 
level of the source lead. The tunneling is due to equilib- 
rium fluctuations between source and QDl. A measure- 
ment of the count rate as a function of /ii provides a way 
to determine both the tunneling rate Ts and the electron 
temperature in the source lead [23 . The situation is re- 
versed at point II in Fig. [TJb). Here, electron tunneling 
occurs between QD2 and the drain, thus giving an inde- 
pendent measurement of Fd and the electron tempera- 
ture of the drain lead. At point III within the triangle 
of Fig. [ijb), the levels of both QDl and QD2 are within 
the bias window and the tunneling is due to sequential 
transport of electrons from the source lead into QDl, 
over to QD2 and finally out to the drain. The electron 



flow is unidirectional and the count rate relates directly 
to the current flowing through the system [17^. Between 
the triangles, there are broad, band-shaped regions with 
low but non-zero count rates where sequential transport 
is expected to be suppressed due to Coulomb blockade 
[cases IV and V in Fig. [TJb,c)]. The finite count rate in 
this region is attributed to electron tunneling involving 
virtual processes. These features will be investigated in 
more detail in the forthcoming sections. 

To begin with, we use the time-resolved charge de- 
tection methods to characterize the system. Typical 
time traces of the QPC current for DQD configurations 
marked by I and II in Fig. [ijb) are shown in Fig. [2|a). 
The QPC current switches between two levels, corre- 
sponding to electrons entering or leaving QDl (case I) or 
QD2 (case II). The change A/qpc as one electron enters 
the DQD is larger for charge fluctuations in QD2 than 
in QDl. This reflects the stronger coupling between the 
QPC and QD2 due to the geometry of the device. A 
measurement of A/qpc thus gives information about the 
charge localization in the DQD. 

In Fig. [2|b) we investigate the charge localization in 
more detail by plotting the absolute change in QPC cur- 
rent A/qpc for the same set of data as in Fig.jlja). The 
detector essentially only measures two different values of 
A/qpc; either A/qpc ~ — 0.3nA or A/qpc ~ — 0.6nA. 
Comparing the results of Fig. |2|b) with the sketches in 
Fig. flic), we see that regions with high A/qpc match 
with the regions where we expect the counts to be due 
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FIG. 2: (a) Typical time traces of the QPC current from con- 
figurations I and II in Fig. [T]^b). (b) Change of QPC current 
A/qpc as one electron enters the DQD, extracted from the 
same set of data as shown in Fig. [l] The two levels corre- 
spond to the QPC detector registering electron tunneling in 
QDl and QD2, respectively. The energy level diagrams de- 
scribe the hole and the electron cycle of sequential transport 
within the finite bias triangles. 
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to electron tunneling in QD2, while the regions with low 
A/qpc come from electron tunneling in QDl. 

The regions inside the bias triangles are described in 
detail in the energy level diagrams of Fig.|2];b). We as- 
sume each QD to hold n and m electrons, respectively. In 
the lower triangle, the current is carried by a sequential 
electron cycle. Starting from the (n,m)-configuration, an 
electron will tunnel in from the source lead at a rate 
Fs making the transition (n, m) (n + l,m). The 
electron then passes on to QD2 at a rate Fc ~ t/h 
[(n + 1, m) (n, m + 1)] before leaving to drain at the 
rate Fd [(n, m + 1) ^ (n, m)]. Since the rate Fc is much 
faster than the detector bandwidth (and Fc ^ Fs, Fc ^ 
Fd), the detector will only register transitions between 
the two states (n, m) and (n,m + 1). Therefore, we ex- 
pect the step height A/qpc within the lower triangle to 
be equal to A/qpc measured for electron fluctuations in 
QD2, in agreement with the results of Fig. [2| 

For the upper triangle, the DQD holds an additional 
electron and the current is carried by a hole cycle. Start- 
ing with both QDs occupied [(n+l,m+l)], an elec- 
tron in QD2 may leave to the drain [(n + l,m+ 1) 
(n + l,m)], followed by a fast interdot transition from 
QDl to QD2 [(n + l,m) (n,m+l)]. Finally, an 
electron can tunnel into QDl from the source lead 
[(n, m + 1) ^ (n + 1, m + 1)]. In the hole cycle, the de- 
tector is not able to resolve the time the system stays in 
the (n + l,m) state; the measurement will only register 
transitions between (n + 1, m + 1) and (n, m + 1). This 
corresponds to fluctuations of charge in QDl, giving the 
low value of A/qpc in Fig. [2|b). Finally, we note that at 
the transition between regions of low and high A/qpc the 
electron wavefunction delocalizes onto both QDs. This 
provides a method for determining the interdot coupling 
energy t [2?. From the data in Fig. [2|b) we find tunnel 
couplings in the range of 10 — 50 /ieV. 
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FIG. 3: Tunneling rates for electrons entering and leaving the 
DQD, measured while keeping the potential of QD2 fixed and 
sweeping the electrochemical potential of QDl. The data is 
measured in a configuration similar to going along the dashed 
line in Fig.[T]^b). The dotted lines are tunneling rates expected 
from sequential tunneling, while the dashed line is a fit to the 
cotunneling model of Eq. ([T]). The solid line corresponds to 
the model involving molecular states [Eq. Parameters 
are given in the text, (b) Schematic drawings of the DQD en- 
ergy levels for three different configurations in (a). At point 
A, electrons tunnel sequentially through the structure. Mov- 
ing to point B, the energy levels of QDl are shifted and the 
electron in QDl is trapped due to Coulomb blockade. Elec- 
tron transport from source to QD2 is still possible through 
virtual processes, but the rate for electrons entering the DQD 
drops substantially due to the low probability of the virtual 
processes. At point C, the next level of QDl is brought inside 
the bias window and sequential transport is again possible. 



III. COTUNNELING 

We now focus on the regions of weak tunneling oc- 
curing in regions outside the boundaries expected from 
sequential transport. In case IV, the electrochemical po- 
tential of QDl is within the bias window, but the po- 
tential of QD2 is shifted below the Fermi level of the 
source and not available for transport. We attribute the 
non-zero count rate for this configuration to be due to 
electrons cotunneling from QDl to the drain lead. The 
time-energy uncertainty principle still allows electrons to 
tunnel from QDl to drain by means of a higher order pro- 
cess. In case V, the situation is analogous but the roles 
of the two QDs are reversed; electrons cotunnel from the 
source into QD2 and leave sequentially to the drain lead. 

To investigate the phenomenon more carefully, we mea- 
sure the rates for electrons tunneling into and out of the 
DQD in a configuration similar to the configuration along 
the dashed line in Fig. [ijb). The line corresponds to 
keeping the electrochemical potential of QD2 fixed within 



the bias window and sweeping /ii . The data is presented 
in Fig. [3j In the region marked by A in Fig. [Sj elec- 
trons tunnel sequentially from source into QDl, relax 
from QDl down to QD2 and finally tunnel out from QD2 
to the drain lead. Proceeding from region A to region 
B, the electrochemical potential /ii is lowered so that an 
electron eventually gets trapped in QDl. At point B, the 
electrons lack an energy 5a = 1-12 — Mi to leave to QD2. 
Still, electron tunneling is possible by means of a virtual 
process [5 . Due to the energy-time uncertainty princi- 
ple, there is a time-window of length ^?i/5a within which 
tunneling from QDl to QD2 followed by tunneling from 
the source into QDl is possible without violating energy 
conservation. An analogous process is possible involving 
the next unoccupied state of QDl, occuring on timescales 
^h/Sb^ where = Eqi — Sa and Eci is the charging en- 
ergy of QDl. The two processes correspond to electron 
cotunneling from the source lead to QD2. Continuing 
from point B to point C, the unoccupied state of QD2 
is shifted into the bias window and electron transport is 
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again sequential. 

In the sequential regime (regions A and C), we fit the 
rate for electrons entering the DQD to a model involving 
only sequential tunneling [dotted lines in Fig. [sj^a)] [24] . 
The fit allows us to determine the tunnel couplings be- 
tween source and the occupied (Fsaj/unoccupied (Fsb) 
states of QD2, giving Fsa = 7.5 kHz, Fsb = 3.3 kHz and 
T = 100 mK. Going towards region B, the rates due to 
sequential tunneling are expected to drop exponentially 
as the energy difference between the levels in QDl and 
QD2 is increased. In the measurement, the rate Fin ini- 
tially decreases with detuning, but the decrease is slower 
than exponential and flattens out as the detuning gets 
larger. This is in strong disagreement with the behavior 
expected for sequential tunneling. Instead, in a region 
around point B we attribute the measured rate Fin to be 
due to electrons cotunneling from source to QD2. 

The rate for cotunneling from source to QD2 is given 
as ^25j : 



t| + Tsb t| + COS (/) VrSa Tgb . . • 

Ob 



(1) 



Here, ta, h are the tunnel couplings between the occu- 
pied/unoccupied states in QDl and the state in QD2. 
The first term describes cotunneling involving the oc- 
cupied state of QDl, the second term describes the co- 
tunneling over the unoccupied state and the third term 
accounts for possible interference between the two. The 
phase (/) defines the phase difference between the two pro- 
cesses. To determine <p one needs to be able to tune 
the phases experimentally, which is not possible from the 
measurement shown in Fig. [sja). In the following we 
therefore assume the two processes to be independent 

= 7r/2). Interference effects between cotunneling pro- 
cesses have been studied in detail in Ref. [15 . 

The dashed line in Fig. [sja) shows the results of 
Eq. ([T]), with fitting parameters = 15 /ieV and = 
33 /ieV. These values are in good agreement with values 
obtained from charge localization measurements. The 
values for Fsa and Fsb are taken from measurements in 
the sequential regimes. We emphasize that Eq. ([T]) is 
valid only if , ^6 ^ , h and if sequential transport is 
sufficiently suppressed. The data points used in the fit- 
ting procedure are marked by filled squares in the figure. 
It should be noted that the sequential tunneling in re- 
gion C prevents investigation of the cotunneling rate at 
small ^5. This can easily be overcome by inverting the 
DQD bias. The rate for electrons tunneling out of the 
DQD [F 

out in Fig. [3]^ a)] shows only slight variations over 
the region of interest. This is expected since fi2 stays 
constant over the sweep. The slight decay of Fout with 
increased detuning comes from tuning of the tunnel bar- 
rier between QD2 and the drain [26] . 

The cotunneling may be modified by the existence of 
a near-by QPC. If the QPC were able to detect the pres- 
ence electron in QD2 during the cotunneling we would 
expect this to influence the cotunneling process. For the 



measurements in Fig.jsja) the QPC current was kept be- 
low 10 nA. This gives an average time delay between two 
electrons passing the QPC of e/igpc ~16 ps. Since this 
is larger than the typical cotunneling time, it is unlikely 
that the electrons in the QPC are capable of detecting 
the cotunneling process. The influence of the QPC may 
become important for larger QPC currents. However, 
when the QPC bias voltage is larger than the detun- 
ing (e^QPc > the fluctuations in the QPC current 
may start to drive inelastic charge transitions between 
the QDs [151 EO] • Such transitions will compete with the 
cotunneling. For this reason it was not possible to ex- 
tract what effect the presence of the QPC may have on 
the cotunneling process. 



IV. MOLECULAR STATES 

The overall good agreement between Eq. ([T]) and the 
measured data demonstrates that time-resolved charge 
detection techniques provide a direct way of quantita- 
tively using the time-energy uncertainty principle. How- 
ever, a difficulty arises as (5 ^ 0; the cotunneling rate 
in Eq. (IT]) diverges, as visualized for the dashed line in 
Fig. (sja)^^ The problem with Eq. ([T]) is that it only takes 
second-order tunneling processes into account. For small 
detuning S the cotunneling described in Eq. ([l]) must be 
extended to include higher order processes [27 . 

A different approach is to assume the coupling between 
the QDs to be fully coherent and describe the DQD in 
terms of the bonding and antibonding molecular states 
[28] [29) . Both the sequential tunneling and the cotunnel- 
ing can then be treated as first-order tunneling processes 
into the molecular states; what we in Fig. [3] referred to 
as cotunneling would be tunneling into an antibonding 
state. The model is sketched in Fig. [4]^a). The bond- 
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FIG. 4: (a) Cotunneling described using molecular states. 
Due to the large detuning the empty antibonding state is 
mainly localized on QD2, but a small part of the wavefunction 
is still present in QDl which allows an electron to enter from 
the source, (b) The rate for electrons tunneling into the DQD 
(Fin) as a function of DQD detuning da- The figure shows 
the same data as in Fig. [S] but plotted on a log-log scale to 
enhance the features at small detuning. The dashed line is 
the results of the cotunneling model in Eq. ([T]), the solid line 
shows the result of the molecular-state model [Eq. ([5|]. 
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ing state is occupied and in Coulomb blockade. Still, 
an electron may tunnel from drain into the antibonding 
state. Due to the large detuning, the antibonding state 
is mainly located on QD2, the overlap with the electrons 
in the source lead is small and the tunneling is weak. 
Changing the detuning will have the effect of changing 
the shape of the molecular states and shift their weights 
between the two QDs. 

To calculate the rate for electron tunneling from source 
into the molecular state of the DQD as visualized in 
Fig. |4ja), we consider the DQD as a tunnel-coupled two- 
level system containing one electron, isolated from the 
environment. We introduce the basis states {^i, ^2} de- 
scribing the electron sitting on the left or the right QD, 
respectively. The two states are tunnel coupled with cou- 
pling t and separated in energy by the detuning S. The 
Hamiltonian of the system is 



H = 



-S/2 t 
t 6/2 



(2) 



The eigenvectors of the Hamiltonian in Eq. ([2| form the 
bonding and antibonding states of the system. 
The eigenvalues give the energies Eb^ Ea of the two 
states, with 

E^ = -^^/W^, Ea = \V^^+5^. (3) 

Note that at zero detuning there is still a finite level 
separation set by the tunnel coupling. The occupation 
probabilities pb and pa of the two states are determined 
by detailed balance. 



1 



PA 



1 



3 fcn T 



(4) 



To calculate the rate for electrons tunneling from source 
into the antibonding molecular state of the DQD as vi- 
sualized in Fig. [4]^a), we project the thermal population 
Pb, Pa of the molecular states and onto the un- 
perturbed state of QDl, ^1. This gives the probability 
pi of finding an electron in QDl if making a projective 
measurement in the ^1 -basis. The measured rate Fin 
is equal to the probability of finding QDl being empty 
(1 — pi) multiplied with Fs, the tunneling rate between 
the source and the unperturbed state in QDl. 



rs(i 



Pi) = rs (1-(pb*b+Pa*a)-*i) 
aanh(^) 



(5) 



For large detuning, the bonding and antibonding states 
are well localized in QDl and QD2, respectively. Here, 
we should recover the results for the cotunneling rate 
obtained for the second-order process [Eq. ([T])]. First, we 
assume low temperature kBT <C so that the electron 



only populates the bonding ground state {pb = 1 and 
PA = 0): 



Tin = Ts - ( 1 



(6) 



In the limit S t the relation reduces to Fin ^ Ts t'^/S'^ 
and the rate approaches the result of the second-order 
cotunneling processes in Eq. ([T]). The advantage of the 
molecular-state model is that it is valid for any detuning, 
both in the sequential and in the cotunneling regime. 

The solid line in Fig.jsja) shows the results of Eq. ([5|. 
The equation has been evaluated twice, once for the oc- 
cupied [(n,m)] and once for the unoccupied state in QD2 
[(n,m+l)]; the curve in Fig. [sja) is the sum of the two 
rates. The same parameters were used as for the cotun- 
neling fit of Eq. ([1]). The model shows very good agree- 
ment with data over the full range of the measurement. 
To compare the results of the molecular-state and the 
cotunneling model in the regime of small detuning, we 
plot the data in Fig. 3]^ a) on a log-log scale [Fig. |4jb)]. 
For large detuning, the tunneling rate follows the 
predicted by both the molecular-state and the cotunnel- 
ing model. For small detuning, the deviations become 
apparent as the cotunneling model diverges whereas the 
molecular-state model still reproduces the data well. 



V. EXCITED STATES 

So far, we have only considered cotunneling involving 
the ground states of the two QDs. The situation is more 
complex if we include excited states in the model; the 
measured rate may come from a combination of cotun- 
neling processes involving different QD states. To inves- 
tigate the infiuence of excited states experimentally, we 
start by extracting the DQD excitation spectrum using 
finite bias spectroscopy [1 . If the coupling between the 
QDs is weak {tc < AEi, AE2, with AEi^2 being the 
mean level spacing in each QD), the DQD spectrum es- 
sentially consists of the combined excitation spectrum of 
the individual QDs. For a more strongly coupled DQD 
the QD states residing in different dots will hybridize and 
delocalize over both QDs. In this section we consider a 
relatively weakly coupled configuration {t ~ 25 /ieV) and 
assume the excited states to be predominantly located 
within the individual QDs. 

Figure [5] shows a magnification of two triangles from 
Fig. [l][b), measured with both negative and positive 
bias applied across the DQD. Excited states are visible 
within the triangles, especially for the case of positive 
bias [marked with arrows in Fig. [sja)]. Transitions be- 
tween excited states occur along parallel lines at which 
the potential of QDl is held constant; this indicates that 
the excited states are located in QDl. To investigate the 
states more carefully, we measure the separate tunnel- 
ing rates Fin and Font along the dashed lines in Fig. [5j 
The results are presented in Fig. [6j together with a few 
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FIG. 5: Finite-bias spectroscopy of the DQD, taken with 
positive (a) and negative (b) bias. The figures are con- 
structed by counting electrons entering and leaving the DQD. 
Excited states are visible, especially for the positive bias 
data [marked with arrows in (a)]. The data was taken with 
Vdqd-sd = ±500 /iV, Vqpc-sd = 250 /iV. 



sketches depicting the energy level configuration of the 
system. 

We begin with the results for the positive bias case, 
which are plotted in Fig. [6|a). Going along the dashed 
line in Fig.jsja) corresponds to keeping the detuning S be- 
tween the QDs fixed and shifting the total DQD energy. 
The measurements were performed with a small detun- 
ing {S ^ 100 fieV) to ensure that the electron transport 
is unidirectional. Because of this, the outermost parts of 
the traces in Fig. [6]^a) correspond to regions where trans- 
port is due to cotunneling [compare the dashed line with 
the position of the triangle in Fig.jsja)]; the regions where 
transport is sequential are shaded gray in Fig. |6|a). 

Starting in the regime marked by I in Fig. |6ja,c), elec- 
trons may tunnel from source into the ground state of 
QDl, relax down to QD2 and tunnel out to the drain 
lead. Assuming the relaxation process to be much faster 
than the other processes, the measured rates Fin and Fout 
are related to the tunnel couplings of the source and drain 
Tin ^ Ts and Fout = Td- Going to higher gate voltages 
lowers the overall energy of both QDs. At the position 
marked by an arrow in Fig.[6]^a), there is a sharp increase 
in the rate for tunneling into the DQD. We attribute this 
to the existence of an excited state in QDl; as shown in 
case II in Fig. [oj^c), the electron tunneling from source 
into QDl may enter either into the ground (n + 1, m) or 
the excited state (n + l*,m), giving an increase in Fin- 
When further lowering the DQD energy another excited 
state comes into the bias window and Fin increases even 
more [second arrow in Fig. [6]^a)]. The rate for tunneling 
out of the DQD shows only minor variations within the 
region of interest. This supports the assumption that the 
excited states quickly relax and that the electron tunnels 
out of the DQD from the ground state of QD2 

Finally, continuing to the edge of the shaded region 
(Vqi ~ — 9.55mV), the potential of QD2 goes below the 
Fermi level of the drain. Here, electrons get trapped in 
QD2 and the tunneling-out rate drops drastically. At 



the same time. Fin increases; when the electron in QD2 
eventually tunnels out, the DQD may be refilled from ei- 
ther the source or the drain lead. The picture described 
above is repeated in the triangle with hole transport 
(-9.25 mV < Vg2 < -8.9 mV). This is expected, since 
the hole transport cycle involves the same QD states as 
in the electron case. An interesting feature is that Fin 
shows essentially the same values in both the electron 
and the hole cycle, while Fout increases by a factor of 
three. The presence of the additional electron in QDl 
apparently affects the tunnel barrier between drain and 
QD2 more than an additional electron in QD2 affects the 
barrier between QDl and source. 

Next, we move over to the case of negative bias 
[Fig. [6]^b)]. Here, the roles of QDl and QD2 are in- 
verted, meaning that electrons enter the DQD into QD2 
and leave from QDl. Following the data and the argu- 
ments presented for the case of positive bias, one would 
expect this configuration to be suitable for detecting ex- 
cited states in QD2. However, looking at the tunneling 
rates within the sequential region of Fig.|6jb), the rate for 
entering QD2 (Tin) stays fairly constant, while the rate 
for tunneling out decreases at the point marked by the 
arrow. Again, we attribute the behavior to the existence 
of an excited state in QDl. 

The situation is described in sketch HI of Fig. ^c). 
The electrochemical potential of QDl is high enough to 
allow the electron in the (n + 1, m)-state to tunnel out to 
the source and leave the DQD in an excited state (n*, m). 
Since the energy difference £^[(n*,m)] — E[{n + l,m)] is 



hole [b] electron 




FIG. 6: (a, b) Tunneling rates for electrons entering and leav- 
ing the DQD, measured along the dashed lines in Fig.jsja, b). 
In (a), we show the results for positive bias across the DQD, 
in (b) the results for negative bias. The shaded areas mark 
the regions where electron transport is sequential, either in 
the electron or the hole transport cycle. The arrows indi- 
cate the positions of excited states. The data was extracted 
from QPC conductance traces of length T = 5s, taken with 
Vqpc-sd = 250 mV. (c) Schematics of the DQD energy con- 
figuration at three different positions in (a, b). 
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FIG. 7: (a) Tunneling rates for electrons entering and leaving 
the DQD, extracted from the same set of data as used in 
Fig.jsjb). The data was measured with Vdqd-sd = —500 /iV. 
The solid lines mark the position of the finite-bias triangles. 
The plot region in the right-hand panel has been extended 
to include the regime investigated in Fig. Is] (b) Energy-level 
diagrams for the two positions marked in (a). In case I, the 
cotunneling itself is elastic, with energy relaxation occurring 
after the cotunneling has taken place. In case II, inelastic 
cotunneling processes are possible. 

smaller than ^[(n, m)] — ^[(n + 1, m)], the transition in- 
volving the excited state appears below the ground state 
transition. As the overall DQD potential is lowered, the 
transition energy involving the excited state goes below 
the Fermi level of the drain, resulting in a drop of Tout as 
only the ground state transition is left available. Similar 
to the single QD case [^, the tunneling-in rate samples 
the excitation spectrum for the (n+l,m)-configuration, 
while the tunneling-out rate reflects the excitation spec- 
trum of the (n, m)-DQD. 

To conclude the results of Fig. [6) we find two ex- 
cited states in QDl in the (n + l,m) configuration with 
AE^ = 180 /ieV and AE^ = 340 /ieV, and one ex- 
cited state in QDl in the (n, m) configuration, with 
A^i = 220 jaeV. No clear excited state is visible in QD2. 
This does not necessarily mean that such states do not 
exist; if they are weakly coupled to the lead they will 
only have a minor influence on the measured tunneling 
rates. Excited states in both QDs have been measured 
in other configurations; there, we find similar spectra of 
excited states for both QDs. 

VI. INELASTIC COTUNNELING 

Next, we investigate the cotunneling process in the 
presence of excited states. Looking carefully at the lower- 
right regions of the negative-bias triangles in Fig. [sjb) , 
we see that the count rates in the cotunneling regions 
outside the triangles are not constant along lines of fixed 
detuning (corresponds to going in a direction parallel to 
the dashed line). Instead, the cotunneling regions seem 
to split into three parallel bands. 

In Fig. [7|^a), we plot the tunneling rates Fin and Fout 
for electrons entering and leaving the DQD, extracted 
from the same set of data as used in Fig. jsjb). The 



thick solid lines mark the edges of the finite-bias triangles. 
Again, the cotunneling rates outside the triangles are not 
uniform; parallel bands appear in Fin for the position 
marked by I and in Fout for the position marked by II in 
the figures. 

To understand the data we draw energy level diagrams 
for the two configurations [see Fig. [7|^b)]. Focusing first 
on case I, we see that the electrochemical potential of 
QDl is within the bias window, whereas QD2 is detuned 
and in Coulomb blockade. The cotunneling occurs via 
QD2 states; electrons cotunnel from drain into QDl, fol- 
lowed by sequential tunneling from QDl to the source 
lead. The picture is in agreement with what is measured 
in Fig.[7|(a); the cotunneling rate (Fin) is low and strongly 
depends on detuning (5, while the sequential rate Fout is 
high and essentially independent of detuning. The three 
bands seen in Fin occur because of the excited states in 
QDl; depending on the average DQD energy, electrons 
may cotunnel from drain into one of the excited states, 
relax to the ground state and then leave to the source 
lead. The state of QD2 remains unaffected by the cotun- 
neling process. For this configuration, we speak of elastic 
cotunneling. 

The situation is different in case II. Here, cotunnel- 
ing occurs in QDl as electrons tunnel directly from QD2 
into the source lead. This means that Fin is sequential 
while Fout describe the cotunneling process. As in case I, 
the cotunneling rate Fout splits up into three bands; we 
attribute this to cotunneling where the state of QDl is 
changed during the process. QDl ends up in one of its 
excited states. The energy of the electron arriving in the 
source lead is correspondingly decreased compared to the 
electrochemical potential of QD2. Here, the cotunneling 
is inelastic. 

The inelastic cotunneling is described in greater detail 
in Fig. [sj In Fig. (sja) we plot the count rate for positive 
and negative DQD bias, measured along the dashed line 
at the right edge of Fig. [TJ^a). Figure [sj^b) shows energy 
level diagrams for negative bias at two positions along 
this line. The bias voltage is applied symmetrically to the 
DQD, meaning that the Fermi levels in source and drain 
leads are shifted by ±6^/2 relative to the Fermi energy at 
zero bias [dotted line in Fig.jsjb)]. In the measurement of 
Fig.jsja) we sweep the average DQD energy while keeping 
the detuning S constant. The average DQD energy is 
defined to be zero when /i2 aligns with the zero-bias Fermi 
energy in the leads [i.e. when /i2 = (/^s + Md)/2]. 

Starting in the configuration marked by A, cotunnel- 
ing is only possible involving the QD2 ground state. Co- 
tunneling is weak, with count rates being well below 1 
count/s. Continuing to case B, we raise the average DQD 
energy. When the electrochemical potential of QD2 is 
sufficiently increased compared to the Fermi level of the 
source, inelastic cotunneling becomes possible leading to 
a sharp increase in count rate. The process is sketched 
in Fig. [sj^b); it involves the simultaneous tunneling of 
an electron from QD2 to the first excited state of QDl 
with an electron in the QDl ground state leaving to the 
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source. The process is only possible if 



S-AE'^=^i-fi2- A£f > /xs - Ml- 



(7) 



Here, A£^f is the energy of the first excited state in QDl. 
The position of the step in Fig. [sj^a) directly gives the 
energy of the first excited state, and we find AE^ = 
180 /ieV. 

Further increasing the average DQD energy makes an 
inelastic process involving the second excited state in 
QD2 possible, giving AE^ = 340 /ieV. Finally, as the 
DQD energy is raised to become equal to half the ap- 
plied bias, the electrochemical potential of QD2 aligns 
with Fermi level of the drain lead. Here electron tun- 
neling mainly occurs due to equilibrium fluctuations be- 
tween drain and QD2, giving a sharp peak in the count 
rate. The excited states energies extracted from the in- 
elastic cotunneling give the same values as obtained from 
finite-bias spectroscopy within the triangles, as described 
in the previous section. The good agreement between the 
two measurements demonstrates the consistency of the 
model. 

The dashed line in Fig. [sja) shows data taken with re- 
versed DQD bias; for this configuration the Fermi levels 
of the source and drain leads are inverted, the electrons 
cotunnel from source to QD2 and the peak due to equi- 
librium tunneling occurs at /i2 = = —300 /ieV. 
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FIG. 8: (a) Electron count rate along the dashed line in 
Fig.[7|^a), measured for both positive and negative DQD bias. 
In the trace, the detuning 6 stays constant and we sweep the 
average DQD energy. The DQD energy is defined from the 
position where the electrochemical potential of QD2 is right 
in the middle between the Fermi levels of the source and drain 
leads [see the dotted line in the energy level diagram in (b)]. 
The steps in the count rate are due to the onset of inelastic 
cotunneling processes in QDl. The data was extracted from 
traces of length T = 10 s, measured with Vqpc-sd = 200 /xV. 
(b) Energy level diagrams for the two configurations marked 
in (a). 
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FIG. 9: Fano factor for electron transport in the cotunnel- 
ing regime of Fig. |3] The data was extracted from traces of 
length T = 30 s. The solid line is the result of Eq. ([8|, which 
assumes independent tunneling events. The minima in Fano 
factor occur at positions where the tunneling rates Fin and 
Tout are equal (see Fig.jsJ. The error bars show standard er- 
ror, extracted by splitting the data into six subsets of length 
T = 5s and evaluating the noise for each subset. 



VII. NOISE IN THE COTUNNELING REGIME 

Using time-resolved charge detection methods, we can 
extract the noise of electron transport in the cotunneling 
regime. For a weakly-coupled single QD in the regime of 
sequential tunneling, transport in most configurations is 
well-described by independent tunneling events for elec- 
trons entering and leaving the QD [TF. The Fano factor 
becomes a function of the tunneling rates [30 : 



Fo = 



F^ 

out 



2e/ (Fi, 



= -(1 
2 



(8) 



with a = (Fin— Fout)/(rin+rout)- For symmetric barriers 
(a = 0), the Fano factor is reduced to 0.5 because of an 
increase in electron correlation due to Coulomb block- 
ade. In the case of cotunneling, the situation is more 
complex. As described in the previous section, cotunnel- 
ing may involve processes leaving either QD in an excited 
state. The excited state has a finite lifetime Trei; during 
this time, the tunneling rates may be different compared 
to the ground-state configuration [8|. We therefore ex- 
pect that the existence of an electron in an excited state 
may induce temporal correlations on time scales on the 
order of r^ei between subsequent cotunneling events. In 
this way, the noise of the cotunneling current has been 
proposed as a tool to probe excited states and relaxation 
processes in QDs [3TJ |32] . 

In Fig. |9j we plot the Fano factor measured from the 
same region as that of Fig. [3j The Fano factor was ex- 
tracted by measuring the distribution function for trans- 
mitted charge through the system [16 . The solid line 
shows the result of Eq. ([s]), with tunneling rates extracted 
from the measured traces. In the outermost regions of 
the graph, the electrons tunnel sequentially through the 
DQD. Here, the Fano factor is reduced due to Coulomb 
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blockade, similar to the single QD case. At the edges 
of the cotunneling regions, the Fano factor drops fur- 
ther down to F = 0.5. This is because the injection rate 
Fin drops drastically as sequential transport becomes un- 
available, while Fout stays approximately constant. At 
some point we get Fin = Tout, which means that the 
asymmetry a is zero and the Fano factor of Eq. (|8| shows 
a minimum. Further into the cotunneling region, the 
Fano factor approaches one as transport essentially be- 
comes limited by a single rate; the cotunneling rate (Fin) 
is two orders of magnitude smaller than the sequential 
rate Font- 

We do not see any major deviation from the results of 
Eq. ([8|, which is only valid assuming independent tun- 
neling events. We have performed similar measurements 
in several inelastic and elastic cotunneling regimes, with- 
out detecting any clear deviations from Eq. ([8|. As it 
turns out, there are two effects that make it hard to de- 
tect correlations due to the internal QD relaxations. For 
the first, the correlation time is essentially set by the re- 
laxation time Trei, which typically occurs on a ~ 10 ns 
timescale. This is several orders of magnitude smaller 
than a typical tunneling time of ~ 1/Fin ~ 100ms j35] . 
Secondly, the slow cotunneling rate limits the amount of 
experimental data available within a reasonable measure- 



ment time. This explains the large spread between the 
data points in Fig.|9]in the cotunneling regime. We con- 
clude that the measurement bandwidth currently limits 
the possibility of examining correlations in the cotunnel- 
ing regime using time-resolved detection techniques. A 
higher-bandwidth detector would solve both the above 
mentioned problems. It would allow a general increase 
of the tunneling rates in the system, which would both 
decrease the difference between Tcot and r^ei as well as 
provide faster acquisition of sufficient statistics. 

To conclude, we have used time-resolved charge- 
detection techniques to investigate tunneling of single 
electrons involving virtual processes. The measurement 
method provides precise determination of all coupling en- 
ergies, which allows a direct comparison with tunneling 
rates expected from time-energy uncertainty. The results 
give experimental confirmation of the equivalence be- 
tween cotunneling through atomic states and sequential 
tunneling into molecular states. In the high-bias regime, 
we measure inelastic cotunneling due to virtual processes 
involving excited states of the double quantum dot. For 
future experiments with a high-bandwidth detector, the 
method may provide a way to probe relaxation processes 
and internal charge dynamics in quantum dots. 
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